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Abstract 

We show how a nonlinear chaotic system, the parametrically kicked nonlinear 
oscillator, may be realised in the dynamics of a trapped, laser-cooled ion, 
interacting with a sequence of standing wave pulses. Unlike the original optical 
scheme [G.J. Milburn and C.A.Holmes, Phys. Rev A, 44, 4704, (1991) ], the 
trapped ion enables strongly quantum dynamics with minimal dissipation. 
This should permit an experimental test of one of the quantum signatures 
of chaos; irregular collapse and revival dynamics of the average vibrational 
energy. 
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It is now well understood that the quantum dynamics of a classically chaotic system will 
show major departures from the classical motion on a suitable time scale [0,3 . Systems 
which classically show chaotic diffusion of a slow momentum-like variable, will cease to 
show diffusion after the 'break-time'. In this region, the momentum distribution is found 
to be exponentially localised, rather than Gaussian. One system which demonstrates this 
phenomenon is a laser cooled atom moving under the dipole force of a standing wave with 
periodically modulated nodal positions [§,§J Recently this prediction achieved experimental 
confirmation in the work of Raizen and coworkers || . The quantum and classical dynamics 
can also differ through the existence of quantum tunnelling between fixed points of the 
classical Poincare section, a phenomenon which has not yet been observed. In this paper 
we consider yet another way in which quantum and classical motion can differ for systems 
which classically do not show diffusion. In this case the transition to chaos in the classical 
system is manifested in the quantum case by a transition from regular collapse and revivals 
of average values with time, to a highly irregular collapse and revival sequence ||. An 
example of such a system is the parametrically kicked nonlinear oscillator |7[. A nonlinear 
oscillator, in which the frequency is linearly dependent on energy, is periodically kicked by 
a parametric amplification process, which momentarily turns the origin into a hyperbolic 
fixed point. The classical dynamics of this system exhibits regions of regular and chaotic 
motion in the Poincare section. 

In reference it was proposed to realise this system using a combination of a Kerr 
optical nonlinearity and an optical parametric amplifier. Unfortunately the required optical 
nonlinearity is usually rather small and accompanied by a large amount of dissipation. For 
this reason, it is unlikely that an optical version will ever show interesting quantum fea- 
tures. However the system does show some interesting classical chaotic behaviour including 
a strange attractor in which unstable periodic orbits can be stabilised by OGY control 
In this paper we show that a much more promising realisation of the quantum version of 
this system may be achieved in the dynamics of a single trapped laser cooled ion interacting 
with a sequence of pulsed optical fields. Furthermore the dynamics may be monitored with 
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high quantum efficiency using the QND vibrational energy measurement scheme proposed 
by Filho and Vogel together with the now standard method of detection on a strong 
probe transition. 

We consider the vibrational motion of a single trapped ion. In the absence of laser 
pulses, the ion executes simple harmonic motion in the quadratic trap potential. As we 
show below, a suitable sequence of laser pulses enables us to realise the nonlinear quantum 
dynamics specified by the map 

\lp n+ l) = U NL U PA \lpn) (1) 

where 

U NL = exp (-iOa^a - i|(a f ) 2 a 2 ^) (2) 

describes an intensity dependent rotation in the oscillator phase plane while 

f/ PA = expg((at) 2 -a 2 )) (3) 

describes a parametric amplification kick. In these expressions the operators a, are the 
ladder operators for the quantised harmonic motion of the trapped ion, which may be written 
in terms of the dimensionless position and momentum operators Xi, X2, by a = X\ + 1X2- 



I. CLASSICAL DYNAMICS 

The classical model corresponding to Eqfll]) is closely related to the model in reference 

0, 

X[ = g (cos(# + ^iR 2 )X l + sin(# + ^R 2 )X 2 ) (4) 
X' 2 = - (cos(# + fiR 2 )X 2 - sin(fl + ^R 2 )X X ) (5) 

where R 2 = X 2 + X 2 and the parametric gain is g = e r . The parameter /1 is simply a scaling 
parameter and the crucial control parameters are the parametric gain g and 9. When 9 = 2tt, 
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easily achievable in experiment as we show below, this model reduces to that discussed in 
reference 0. The origin is then a saddle. If however Icoshrcos^l < 1, where r = Ing, 
the origin is elliptically stable. The bifurcation diagram is shown in figure [I]. The origin is 
elliptically stable in the regions bounded by the solid curves 9 = iarccos ( eos 1 h - j + nir and 
labelled E * *. Away from the origin an infinite number of period- 1 orbits lie on lines with 
a slope of ±e(~ r \ independent of 9 as in the original model. However the radius at which 
the period-1 orbits lie does depend on 9, and is determined by the equation 

- tan(0 + fiR 2 ) = ± sinh r (6) 

As 9 is increased the period-1 orbits move inwards collapsing on the origin and changing its 
stability. Those in the odd quadrants are always unstable. In the even quadrants a finite 
number (maybe zero) lying closest to the origin are stable depending on whether 

1 > fxR 2 tan(0 + nR 2 ) > (7) 

They loose stability on the curves 9 = arctan(sinhr) — -r4- + nn, shown as dashed (for 
n=-l,0,l) in figure [1]. The various parameter regions are labelled with a sequence of Es and 
Hs standing for elliptic or hyperbolic. The first letter refers to the stability of the origin 
and the subsequent letters to the stability of the pairs of period-1 orbits in the 2nd and 
4th quadrants starting with those closest to the origin and moving out. Since the elliptic 
period-1 orbits always occur closest to the origin in the region EEH,for instance, only one 
pair of period-1 orbits are elliptically stable. 

In figure |2] we show typical Poincare sections of the classical dynamics for various values 
of g, and 9 = 2tt. Note the transition to almost global chaos for large values of g. 

II. QUANTUM DYNAMICS 

To compare the quantum and classical dynamics we compute the mean energy as a 
function of the number of kicks. In the quantum case we assume the oscillator is initially 



4 



prepared in the vibrational ground state, while in the classical case we use an initial Gaussian 
phase-space distribution with the same moments as the initial quantum state. The detailed 
results are given in reference |7], which we summarise here. In the classical case, the average 
energy saturates at an almost time independent value, as the initial distribution becomes 
spread around the period-one fixed points. The time taken to reach the saturation value 
is longer for an initial state in the regular region, than for one located in a chaotic region. 
In the quantum case, the mean energy follows the classical for a short time, but instead of 
saturating at an almost stationary value, continues to oscillate with a characteristic collapse 
and revival envelope. The transition from regular to chaotic dynamics in the classical system 
is marked in the quantum system by a transition from a regular collapse and revival sequence 
to a highly irregular and more rapid collapse and revival sequence. This is due to the fact that 
a state localised in the chaotic region has support on very many quasi-energy eigenstates, 
compared to a state localised in a regular region. 

We now show how this system may be realised in a trapped ion experiment. Following 
Filho and Vogel || we first consider an ion trapped at an antinode of an optical standing 
wave tuned to the atomic frequency, the carrier frequency. In an interaction picture at 
frequency v the interaction Hamiltonian is, 

Hi = -Mlifa}ao x + h—t-(a v ) 2 a 2 <j x (8) 

where Q is the Rabi frequency and 77 is the Lamb-Dicke parameter. If the ion is first prepared 
in an eigenstate of a x by the application of a 7r/2 pulse, it will remain in this state. In that 
case the vibrational motion experiences a linear and a nonlinear phase shift. The interaction 
is precisely the kind required to realise Unl, with 6 = TVtrj 2 /i = ^y 2 - where T is the time 
for which the standing wave pulse is applied. 

In order to realise the parametric kick we use a Raman interaction with two laser pulses, 
one tuned to the first lower sideband, the other to the first upper sideband. This scheme 
has been successfully used by Wineland and coworkers to prepare the ion in a squeezed 
vibrational state |IJ. In this case the interaction Hamiltonian, in an interaction picture at 
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frequency v, is 

H PA = iTik(cl 2 - (a f ) 2 ) (9) 

where k = ^g^ 2 ^ and Qi is the Rabi frequency for each of the Raman pulses, each of 
which is detuned from the atomic transition by Si such that 8\ — 62 = 2u. The Lamb-Dicke 
parameter for the Raman transition is 77^ = 5kxo, where 8k is the wave vector difference 
for the two Raman beams and xq is the rms fluctuations of position in the ground state of 
the trap. Note that this Hamiltonian is independent of the ion electronic state. The ion 
evolves according to the parametric kick transformation with g = e KT for a Raman pulse of 
duration T. In summary, we first prepare the ion electronic state by the application of a 
7r/2 pulse, then apply the pulse at the carrier frequency, followed by the Raman parametric 
kick pulses. Pulses at the carrier frequency and the Raman pulses then alternate to effect 
each kick. 

A particular advantage of this realisation is that it is easy to read out the energy of the 
vibrational motion at any stage. If we reprepare the ion in the ground electronic state at 
the end of the carrier frequency pulse, its evolution will be dominated by the first term in 
Eq(^). This is a QND measurement interaction for the vibrational quantum number. It 
causes the Bloch vector describing the two level system to precess around the x-axis by an 
angle proportional to a^a. The probability to find the atom in the ground state a time r 
after the last squeezing pulse is 

00 

P g (r) = P(n) cos 2 0rn (10) 

n=0 

where = Qi] 2 /2. A set of readouts of the atomic state a time after an appropriate number 
of kicks thus enables a moment of the photon number to be constructed directly. This is 
similar to the scheme used by Wineland and coworkers to observe the Jaynes-Cummings 
collapse and revival sequence UTOf . 

In figures |3] and |] we plot P g {r) = (cos 2 0rn) for the same two initial states and parameter 
values as discussed in reference 0. In figure §, the initial state is localised in a classically 
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regular region of phase space. A regular collapse and revival sequence is evident in the 
quantum case (fig. ^a). In figure || however, we have the case of an initial state localised 
in the chaotic part of the classical phase space, leading to an irregular revival sequence. In 
both cases the corresponding classical moment is shown as well. 

We can estimate reasonable values of the parameters 0,g,fi from recent experiments on 
trapped ions. From Meekhof et al. ||10|| , typical values are Q = 10 6 , rj = 0.2, g = 6. For 
an interaction time of T = 10/xs, we have 9 « 1.0, /i ~ 0.02. It is then apparent that the 
parameter values of interest are achievable in an experiment. 
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FIGURES 

FIG. 1. Bifurcation diagram for the origin and the period-1 orbits close to the origin. The 
bifurcation curves for the origin are represented as solid lines while those for the period-1 orbits 
are dashed. The regions are labeled by a finite sequence of E and Hs, standing for elliptic and 
hyperbolic. The first letter refers to the stability of the origin, the second to that of the pair of 
period-1 orbits closest to the origin in the 2nd and 4th quadrants (those in the 1st and 3rd are 
always hyperbolic), the third to the pair of period-1 orbits next closest to the origin etc. 

FIG. 2. Phase space portraits for the classical map. In all cases [i = 0.017T, 9 = 2tt. (a) g = 1.0, 
(b)g = 1.2, (c) g = 1.5, and (d)g = 2.0. 

FIG. 3. Initial states centered in a classically regular region: (a) Plot of the quantum average 
P 9 (r) = (cos 2 4>Tn) for an initial coherent state, (b) Plot of cos 2 (4>tE) (where E is the classical 
mean energy) for an initial Gaussian distribution. In both cases the initial state is centered at 
(0,0), with g = 1.2, = O.Olvr, 9 = 2vr, and (pr = 0.01. 

FIG. 4. Same as Fig. (|3|) except initial states are centered in a classically chaotic region, at 
(1,0), with g = 1.5. 
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